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Abstract. Using the d^2_y2-\-dz2 two orbital model of the high Tc cuprates obtained from 
the first-principle calculation, we show that the material dependence of the Fermi surface shape 
can be understood by the degree of the mixture between the (1^2 _y2 and the d^2 orbitals. 
We explain, through investigating the tightbinding hopping integrals, why some cuprates have 
square shaped Fermi surface, while others have more rounded ones. From this viewpoint, we 
explain the experimentally observed correlation between the curvature of the Fermi surface and 



1. Introduction 

In the high Tc cuprates, the 'main band' having strong Cu3(ia,2_^2 orbital character constructs 
the Fermi surface, and the single band model that considers only the main band has often 
been adopted in the theoretical studies. Such single band or related models have provided 
us many understandings, but there still remain unresolved problems. One of the issues often 
discussed with controversy is the relationship between the curvature of the Fermi surface and the 
critical temperature (Tc). It is well known that even within the single-layered cuprates, there is 
significant difference of Tc, for example, the La compound with Tc ^ 40K and the Hg compound 
with Tc ^ 90K. From the band structure point of view, the La cuprate has relatively square 
shaped (diamond-like) Fermi surface, while the Hg material has a more round one. In fact, it 
has been recognized that low Tc materials have square shaped Fermi surface, while high Tc ones 
tend to have round Fermi surface [U |2]. Although some phenomenological[3j or t — J model [HE] 
studies describe such tendency, a number of Hubbard-type-model studies with realistic values of 
the on-site U have not succeeded in reproducing such a tendency [6J. For example, the dynamical 
cluster approximation studies that adopt a single band model [7J , or a more realistic three band 
model [8J that considers Pcr orbitals show the opposite tendency. 

To give insight into this long standing problem, we have constructed a dx2_y2-\-d^2 two orbital 
model that considers all the orbitals that have eg symmetry [9J. In this proceedings, we focus on 
the relationship between the conventional single orbital model and the two-orbital model, and 
explain the difference of the parameters that gives the material dependence. 



2. d^2_y2 + d^2 two-orbital model 

In Fig.l, we show the band dispersion of the two-orbital model constructed by exploiting 
maximally localized Wannier orbitals,[10j which are obtained from the first-principles calculation 
results[llj. The thickness of the lines represents the strength of the respective orbital character. 
In the La compound, a significant amount of d^2 character is present in the main band around 
the wave vector k = (tt, 0)(0, tt) (denoted as N in La and M in Hg), while in the Hg compound 
such a mixture of the d^2 component is absent [121 [13 [HI [IS]- 
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Figure 1. The band structure of the two orbital model for La2Cu04 (left) and 
HgBa2Cu04 (right). The insets depict the Fermi surfaces (for the total band filling n — 2.85). 
The thickness of the lines represents the strength of the respective orbital character. 



3. Comparison from the viewpoint of the hopping integrals 

Let us now look into the hopping integrals of the models. In table I, the hopping integrals within 
the dx2_y2 orbitals of the two orbital model are displayed for the La and Hg compounds. For 
comparison, we also show the values of the single band model obtained by the same method. 
Conventionally, the difference of the Fermi surface curvature is represented by the value of the 
second nearest neighbor hopping t2 and the third nearest hopping ts, where large \t2\ and 1^3 1 
gives more rounded Fermi surface. This tendency is in fact seen in the table. 

Table 1. Hopping integrals within the dx2_y2 orbital for the single and two-orbital models, and 
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Two-orbital 






La 


Hg 


La 


Hg 


ti[eV] 


-0.444 


-0.453 


-0.471 


-0.456 


t2[eV] 


0.0284 


0.0874 


0.0932 


0.0993 


t3[eV] 


-0.0357 


-0.0825 


-0.0734 


-0.0897 


(I*2| + |t3|)/|tl| 


0.14 


0.37 


0.35 


0.41 


A^[eV] 






0.91 


2.19 



If we turn to the hopping integrals of the two orbital model, the values shown in Table 
I is surprising in that the drj.2_y2 orbital has large \t2\ and 1^3 1 even for the La compound, 
just as in the Hg compound. Namely, the curvature of Fermi surface is not governed by the 
drj.2_y2 distant hoppings in the two orbital model that considers the d^2 orbital explicitly, and 
another parameter plays an important role. In fact, the parameter is the on-site energy difference 
/S.E = E^2_y2 —E^2. In other words, AE in the two-orbital model determines the (|t2| + |t3|)/|ti| 
ratio in the effective single-orbital model. This can be understood as follows. 
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Figure 2. A schematic image of the (1^2 _y2 -(1^2 _y2 diagonal hopping within the single (a) and 
two-orbital(b) models. The left side of (b) shows the direct path, while the right side is the 
indirect path via the (1^2 orbital. 



An electron can hop from a dj,2_y2 orbital to its second nearest neighbor mainly via the 



direct diagonal drj.2_y2-drj.2_y2 path or the indirect drj.2_ 



d^ 



drj.2_y2 path. The latter can 



be considered as a second order perturbation process when the d^2 orbital degrees of freedom 
is integrated out to obtain the effective single band model. As shown in table I, the value of 
the direct hopping is nearly equal between the La and the Hg compound, so that the material 
dependence comes from the indirect term. 

The amplitude of the indirect term is generally large for smaller /S.E. In the La compound, 
the relatively small AE" mainly enhances the indirect path drj.2_y2 d^2 drj.2_y2. From the 
first-principles result, the direct hopping and the indirect term are found to have the opposite 
sign, so the cancellation is strong when the amplitude of the indirect path is strong. As a result, 
the small for the La compound gives the small effective ^2 and t^. 

The above analysis shows that the curvature of the Fermi surface in the La cuprate is 
suppressed due to the large mixture between drj.2_y2 and d^2 orbitals. Therefore, the strong 
(i^2 mixture in the Fermi surface around the wave vectors (tt, 0)/(0, tt) is the origin of the weak 
curvature. 



4. The relationship between Fermi surface and Tc 

Now, we finally discuss how AE" affects the many body properties, especially rf-wave 
superconductivity. We consider a many body Hamiltonian that considers the on-site multiorbital 
interactions with realistic values. In Fig. 3, we plot the rf-wave eigenvalue of the linearized 
Eliashberg equation, where the Green's functions are obtained within the fluctuation exchange 
approximation [ini [H]. The eigenvalue reaches unity when T — Tc, so A at a fixed temperature 
(here T = 0.01 [eV]) can be used as a qualitative measure for Tc. In the figure, the points indicated 
by the arrows are the results of the model constructed from the experimentally determined lattice 
parameters [T5t [T9] . This shows that the Hg compound has higher Tc than the La compounds, so 
the result is consistent with the experiments. We also plot A when is varied hypothetically in 

d 2 _ 2 d 2 _ 2 

the La system. Note that ~^ and t^"" ~^ are fixed because in the first principles calculation 
they are not so much material dependent, as mentioned above. From FigjSl it can be seen that 
A increases monotonically with AE", and comes close to the value of the Hg result. From this 
result, we can say that AE governs both the shape of the Fermi surface (via the effective ^2 and 
ts) and Tc. 



5. Conclusion 

In conclusion, we have shown that the degree of the orbital mixture controls the shape of the 
Fermi surface, and this gives the correlation between the Fermi surface shape and Tc. In this 
picture, the parameter that largely contributes to the material dependence is A£^, the on- 
site level offset between the (1^2 _y2 and d^2 Wannier orbitals. This picture also explains the 
conventionally adopted material dependence of the second and third nearest neighbor hoppings 
in the single band model. 
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Figure 3. The eigenvalue A of 
the Eliashberg equation for d-wave 
superconductivity plotted against 
AE = E^2_y2 — E^2. The points 
with arrows are the results obtained 
using the original lattice structure 
(determined experimentally). Cir- 
cles are the results obtained by hy- 
pothetically varying AE in the La 
system. 
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